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ABSTRACT 

We use high-resolution hydrodynamic simulations to investigate the density profile 
of hot gas in clusters of galaxies, adopting a variant of cold dark matter cosmologies 
and employing a cosmological N-body/smoothed particle hydrodynamics code to follow 
the evolution of dark matter and gas. In addition to gravitational interactions, gas 
pressure, and shock heating, we include bremsstrahlung cooling in the computation. 
Dynamical time, two-body relaxation time, and cooling time in the simulations are 
examined to demonstrate that the results are free from artificial relaxation effects and 
that the time step is short enough to accurately follow the evolution of the system. 
In the simulation with nominal resolution of 66h~^ kpc the computed cluster appears 
normal, but in a higher (by a factor 2) resolution run, cooling is so efficient that the 
final gas density profile shows a steep rise toward the cluster center that is not observed 
in real clusters. Also, the X-ray luminosity of 7 x lO^^ergss"^ far exceeds that for 
any cluster of the computed temperature. The most reasonable explanation for this 
discrepancy is that there are some physical processes still missing in the simulations 
that actually mitigate the cooling effect and play a crucial role in the thermal and 
dynamical evolution of the gas near the center. Among the promising candidate 
processes are heat conduction and heat input from supernovae. We discuss the extent 
to which these processes can alter the evolution of gas. 

Subject headings: cosmology: miscellaneous — galaxies: clusters: general — 
hydrodynamics — intergalactic medium — methods: numerical 
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1. Introduction 



Hot X-ray-emitting gas in clusters of galaxies contains a variety of information relevant to 
many fields in astrophysics. Density and temperature profiles of the gas give among the most 
reliable estimates of cluster mass, which are of unparalleled importance to cosmology (e.g., Bahcall 



1995 ). They will also reflect physical processes that have played a crucial role in the thermal and 
dynamical evolution of the gas. In addition to gravitation, hydrodynamics, and shock heating, 
such processes may include radiative cooling, heat conduction, and feedback from star formation. 

One of the central questions about the cluster gas structure involves the presence of a core. 
Observations have revealed that the gas density profile has a distinct core, inside of which density 
is nearly constant. Gravitational N-body simulations have demonstrated that gravity alone cannot 
produce such a core in an object formed as a result of hierarchical structure formation; halos 
formed in high-resolution N-body simulations have density profiles with significant slope toward 
the center up to the resolution limit (Navarro, Frenk, and White 1996, 1997; |Fukushige and 



Makino 199^ ; [Moore et al. 1998^ . For example, Navarro, Frenk, and White (1996, 1997) found 
that the density profiles of halos with a wide range of masses can be fitted by 

Pc 



p{r) 



(1) 



{r/rs){l + r/rsr 

(hereafter referred to as the NFW profile), with pc and being the fitting parameters. In the 
NFW profile, gas temperature would approach zero in a cluster's central regions (T oc r) were 
it to have the same profile as dark matter. While this model does produce a convergent X-ray 
luminosity, the temperature structure is different from what is observed. It follows that processes 
other than gravity are responsible for the formation of the cores. Makino, Sasaki and Suto (1998) 
pointed out that the gas distribution develops a core in a dark matter halo having a central cusp 
(e.g., the NFW profile) if the gas is isothermal and in hydrostatic equilibrium. Even if the gas 
obeys a less stringent limit of constant entropy (T jr?!"^ — > const), it would avoid a central cusp 
and have an apparently constant-density, constant-temperature core. But it is not evident what 
physical mechanisms could enforce either isothermality or isentropy. 

Many authors have studied the evolution of clusters using simulations that include 



hydrodynamics ( 


Evrard 199C; ' 


rhomas and Couchman 1992; 


Katz and White 1993 


; Bryan et al. 


19944 1994b; |Kang et al. 199^ 


:; jNavarro, Frenk, and White 1995; 


Bartelmann and Steinmet^ 


1996 


; Bryan and Norman 1997, 1998; Eke, Navarro, and Frenk 1995; Pen 1998; Yoshikawa, Itoh, 


and Suto 1998 


). In general these simulations have succeeded in producing clusters that have 



cores similar to those observed. However, the situation is far from satisfactory for at least two 
reasons. First, artificial two-body relaxation may affect the dynamics, especially in the central 
region. Indeed, Steinmetz and White (1997) showed that this effect gives rise to artificial energy 
transfer from dark matter to gas. Both spatial and mass resolutions are only marginally adequate 
in most published works. Second, almost none of the simulations include radiative cooling of gas. 
Although cooling is probably unimportant in the outer part of a cluster ( Sarazin 1986| ), it may 
affect the dynamics in the central region. 



-3- 



This work is an investigation of cluster gas using hydrodynamic simulations that attempt to 
address these limitations. First, we minimize two-body relaxation effects for a given computational 
cost by employing a multiresolution technique that enables us to improve resolution only inside the 
clusters where we really need high resolution. Second, we include cooling due to bremsstrahlung, 
which dominates cooling of gas at above ^ lO^K. We ignore line cooling, but this is not a serious 
problem because we focus our attention on X-ray-emitting gas. A practical reason for ignoring 
line cooling is that doing so enables us to avoid very short timescales in moderate temperature 
(10'^ - lO^K), high density regions. Allowing for line cooling would strengthen the conclusions of 
this paper. 

We organize the rest of the paper as follows. In § 2 we describe the method and parameters 
of the simulations. We present the results in § 3. As we will see, the most important result is 
that cooling and increased resolution give rise to a density profile of the gas that rises steeply 
toward the center and consequently produces excessive X-ray luminosity. In § 4 we discuss the 
implications of our results in connection with physical processes that are still missing in the 
simulations. In § 5 we give our conclusions. 

2. Simulations 

We work with a cold dark matter-dominated universe with h = 0.6, = 1, A = 0, Ob = 0.07, 
'Aspect = 0.81, and fjg = 0.50, where h is the Hubble constant in units of 100kms~^Mpc~^, f^o 
is the cosmological density parameter, A is the cosmological constant, 0,b is the baryon density 
parameter, ngpect is the power-law index of the primordial density fluctuations, and as is the rms 
density fluctuations on 8^~^Mpc scale. 

These parameters are chosen so that the universe is at least marginally consistent with all of 
the following observations: the amplitude of cosmic microwave background anisotropies observed 
by COBE, the observational bound on the Hubble constant, the abundance of light elements, the 
baryon fraction in clusters, the number density of clusters, and the age of the globular clusters. 
The model is similar to the tilted cold dark matter (TCDM) model of Cen and Ostriker (1993). 

We employ a cosmological N-body -|- smoothed particle hydrodynamics (SPH) code (see e.g. 
Monaghan 1992 for a review of SPH). 

The N-body part of the code uses the Barnes-Hut (1986) tree algorithm and is described in 

Suginohara et al. (1991). The SPH part of the code is described in Suginohara (1992). Overall, 
the code is similar to the one described in Katz, Weinberg, and Hernquist (1996). An application 
of the earlier version of the code is in Suginohara (1994, 1995). 

We carry out a multiresolution simulation that allows us to achieve high resolution only inside 
the clusters. Wc prepare two sets of initial conditions with different mass resolution in a cubic 
box with comoving size L = 42 /i~^Mpc. These sets have exactly the same realization, up to the 
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Nyquist wavenumber corresponding to the lower resolution, of the theoretical power spectrum of 
density fluctuations at z = 19. The number of particles in the higher resolution realization is 
8 times that in the lower resolution one. In each realization we assign to all particles the same 
gravitational softening length e, which is fixed in the comoving coordinates and is equal to 0.1 
times the mean interparticle separation of dark matter particles. The initial gas temperature is 
set to 10^ K. Then we use the following two-step procedure. First, we evolve the system from 
the lower resolution realization. When the low-resolution (LR) run has been carried out up to 
the present epoch we identify cluster centers. Then we come back to the initial distribution of 
particles and identify those regions that will end up with overdensity larger than 140 around the 
cluster centers. We replace the particles in the identified regions with those in the high-resolution 
realization. Then, starting from this initial condition, we run the second simulation, which is 
hereafter referred to as the multiresolution (MR) run. Katz and White (1993) and Navarro, 
Frenk, and White (1996, 1997) used a similar method to construct initial conditions with variable 
resolution. 

The ratio of the number of gas particles to that of dark matter particles is always set as 1/8. 
Thus in either the low-resolution or high-resolution realizations, the masses of gas particles, nig, 
and dark matter particles, m^, are comparable {mg/nid = 0.6). Since two-body relaxation time 
is inversely proportional to max(mg,mrf), it is maximized when mg = rrifi under the condition 
that the available CPU time is fixed. Thus, with the available CPU cost, the artificial two-body 
effect is significantly reduced compared to when we used an equal number of gas and dark matter 
particles, as has been done in most previous work. Indeed, the mass of a dark matter particle 
for the high-resolution realization, 1.5 x lO^'^M©, is less than one-third of the upper bound for 
two-body heating to be subdominant ( Steinmetz and White 1997| ) for the cluster that is mainly 



discussed in the following sections. 

We allow each particle to have its own time step. The time step of each particle satisfies 
At = ( At)sys/2"=*<=p , where (At)sys = 2.4x lO'^yr. The time-step index 

'T'stepj which is a nonnegative 
integer, is chosen so that the following three conditions are satisfied: (i) the Courant condition 
is satisfied with the Courant number taken to be 0.3 ( Katz, Weinberg, and Hernquist 1996D , (ii) 



expected relative change in the internal energy within a time step does not exceed 0.05, and (iii) 
expected change in position within a time step does not exceed 0.05e. The parameters for the 
simulations are summarized in Table 

Bremsstrahlung cooling is included in the simulations. Throughout this paper we assume the 
mass fraction of helium, Y , to be 0.24 and the effect of metals to be negligible, and we use the 
approximation that g = 1.2, where g is the frequency averaged Gaunt factor. Then bremsstrahlung 
emissivity is given by 

Kb = 5.2 X 10~2^^l/2n2ergss-^cm-^ (2) 



where the gas temperature T is in units of K and the gas density n is in units of cm ^ ( Rybicki| 
and Lightman 197S| ). 
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Table 1: Simulation parameters 



Run 


Resolution ^ 


AT/ 


Nf 






e{h -"^kpc)^ 


/imin(/i"^kpc)^ 


LR 


low 


32,768 


262,144 


7.3 X 10^° 


1.2 X 10^^ 


66 


66 


MR 


low 


27,962 


223,696 


7.3 X 10^° 


1.2 X 10^^ 


66 


66 


MR 


high 


38,448 


307,584 


9.2 X 10^ 


1.5 X 10^° 


33 


33 



^In run LR, the entire volume is represented by low-resolution particles. In run MR, those regions 
which will end up in clusters are represented by high-resolution particles, and the remainder of 
the volume by low-resolution particles. 

'-'The number of gas particles. 

"^The number of dark matter particles. 

^The mass of a gas particle. 

'^The mass of a dark matter particle. 

^Gravitational softening length. 

^Minimum smoothing length of gas particles. 



3. Results 

In this section we present results of the simulations. We have obtained 12 clusters with 
M4 > 2 X 10^^ Mq, where Ma is the total mass inside a radius of 1.5h ^Mpc. We have identified 
these clusters in the following way. First, we extract gas particles with gas density larger than 
10~^cm~^ out of the whole simulation. Next we group these gas particles using the conventional 
friends-of-friends algorithm, connecting particles with separation smaller than 0.2 times the mean 
interparticle separation. We find the center of mass of each group. Then we come back to the 
entire simulation and extract a spherical region of a radius 2.5/i~^Mpc, whose center matches the 
center of mass of the identified group. Finally we redefine the cluster center as the position of the 
gas particle that has the largest X-ray emissivity. 

3.1. Density and temperature profiles 

In the remainder of the paper we concentrate, unless otherwise stated, on the richest cluster, 
with Ma = 6.9 X IO^^Mq. The half-mass velocity dispersion is 790kms~^, and the cluster is in 
the range of normal X-ray-emitting clusters. Figure |l| shows the contour of the column density 
distribution of gas and dark matter projected on a plane in runs LR and MR. The difference 
for the two runs is already manifest, especially in the gas distribution, despite the fact that the 
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only difference in tlie two runs is the resolution. In run MR the density profile is much more 
centrally concentrated. Moreover, the cluster in run MR has two distinct subclumps, which are 
not apparent in run LR. 

The difference is also clear in Figure ||, which shows total mass density p, gas density n, 
and temperature T, averaged over a spherical shell, as functions of the radius. In run LR, the 
logarithmic slope of the total mass density profile gradually changes from ~ —3 in the outer part 
to ~ —1 in the inner part, i.e., the simulation fits the NFW profile (eq. 

The gas density profile in run LR does not fit the NFW profile but is well fitted by the 
observationally motivated (3 profile: 



n(r) = no 



-3/3/2 



(3) 



with riQ = 5.4 x 10~'^cm~^, = 140/i~^kpc, and /? = 0.75. The gas is nearly isothermal with 
T = (3-4) X 10^ K in run LR. The X-ray luminosity in this run is 1.9 x lO^^ergss"^ and is in the 
normal range observed. This run has resolution and other properties typical of high-resolution 
numerical simulations of hydrodynamical cosmology. If we had stopped our work at this point, we 
would have concluded that there was a good correspondence between observation and theory. 

In contrast, in the higher resolution run MR the total mass density profile shows a steeper 
rise toward the center than in run LR. The gas density profile also shows a steep rise toward the 
center and cannot be fitted at all by the [3 profile. The temperature in run MR significantly drops 
toward the center. There is a sharp peak in the gas density profile and drop in the temperature 
profile at r ~ 0.3/i~^Mpc; this corresponds to the smaller subclump seen in Figure ||. These results 
in run MR can naturally be interpreted as cooling of gas having decreased the temperature and 
increased the gas density in the central part of the subclumps. 

No observed clusters resemble the simulated cluster in run MR. The X-ray intensity is smooth 
in the central region in roughly half of the observed clusters, which implies that the gas density 
is nearly constant. The remainder of the clusters do show a peak in X-ray emission, which is 
claimed to be a feature of a cooling flow ( [Fabian, Nulsen, and Canizares 1991 ). However, the 



inferred gas density at the center of "cooling-flow" clusters is much smaller than in our simulated 
cluster. There is also a discrepancy between the simulation and the observations in the total X-ray 
luminosity. The bolometric luminosity Lboi of the cluster in run MR is 7.4 x 10^^ergss-\ far too 
high for a cluster with temperature 3keV; typically Lboi ~ 10"^^ ergs s~^ for observed clusters 
with this gas temperature, and none of these have bolometric luminosity exceeding 10^^ ergs s~^ 



QDavid et al. 1993^ . 



We have seen that as a result of increasing the resolution the simulated cluster has moved 
farther from the observed ones. This is probably because cooling becomes too efficient. The origin 
of this discrepancy is discussed in § ^. 

In general, the other 11 clusters in the higher resolution run have steep gas density profiles 
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Fig. 1. — Contours of the column density S of the richest cluster in the simulation, projected 
on a plane. The neighboring contour lines correspond to increments of 0.2 in log E. (a) Column 

density of the gas in run LR. The contour line that is nearest to the cluster center corresponds to 
E = 10~^'^gcm~^. (b) Column density of the dark matter in run LR. The contour line nearest 
to the cluster center corresponds to S = lO""'^''^ gcm~^. (c) Column density of the gas in run 
MR. The contour line nearest to the cluster center corresponds to S = 10~^'^gcm~^. (d) Column 
density of the dark matter in run MR. The contour line nearest to the cluster center corresponds 
to E = 10~^-2gcm-2. 
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0.1 1 

r (h-^Mpc) 

Fig. 2. — Profiles of (a) total mass density, (b) gas density, and (c) temperature in the cluster in 
runs MR {solid line) and LR {dashed line). The abscissa is the physical radius. The solid and dotted 
vertical arrows indicate the gravitational softening length e for runs MR and LR, respectively. 
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in the central part, similar to the richest cluster. The central density of six of them is greater 
than 10~^cm~^, and two of them have bolometric luminosity exceeding lO^^ergss"^ although they 
have velocity dispersion of only 500-600 km s~^, considerably smaller than the richest one. Thus 
the properties of the richest cluster discussed above are not to be attributed to a particular initial 
condition. 



3.2. Evolution 

Next we examine the time evolution of the cluster in run MR. Figure |^ shows the total mass 
density, gas density, and temperature profiles at various redshifts. The temperature increases at 
any radius from z = 1.4 to z = 0.7, but at z ~ 0.7 the temperature in the central part begins to 
decrease, and it keeps falling up to z = 0. The gas density in the central part increases rapidly 
from z ~ 0.7 to z = 0. We also notice that the smaller subclump is approaching the main 
subclump. These results support our view in the previous subsection that cooling has resulted in a 
rise in the gas density profile toward the center. On the other hand, the total mass density profile 
in the central part is already quite steep at z = 1.4. This indicates that the steeper rise in run MR 
than in run LR is a consequence not of the effectiveness of cooling but of the improved accuracy 



in following the evolution of the dark matter. As some authors have pointed out (Fukushige and 



Makino 1997; Moore et al. 199^ ), the "universal" density profile of dark halos claimed by Navarro, 



Frenk, and White (1996, 1997, eq. [1]) is still a controversial issue and may depend on numerical 
resolution. Clearly, our higher resolution run allowed the central density to reach higher levels 
than the lower resolution run even in the absence of cooling. But then, in run MR the density 
reached a level where the cooling time became short enough for the classical isobaric radiative 
cooling instability to begin. Allowance for metal-line cooling would have increased the radiative 
losses by about a factor of 2 and would have accelerated the trends noted. 

We find that the cooling time in run LR is greater than the Hubble time up to the resolution 
limit. Therefore the main reason for the inefficiency of cooling in run LR is simply that the 
achieved density is not high enough, although two-body heating ( Steinmetz and White 1997| ) may 
have played an additional role in preventing the gas from cooling. 



3.3. Timescales 

Now we investigate various timescales in the simulation both as a function of radius and as a 
function of time. The purpose is to confirm that no obvious numerical artifacts are invalidating 
the results of the present simulations and that cooling is indeed what has caused the gas density 
profile to rise steeply toward the center. 

The timescales we are going to discuss are the following: 
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0.1 1 



r (h-^Mpc) 

Fig. 3. — Profiles of (a) total mass density, (b) gas density, and (c) temperature in the cluster in 
run MR at z = 1.4 [dot-dashed line), z = 0.7 {long-dashed line), z = 0.3 {short-dashed line) and 
z = {solid line). The abscissa is the physical radius. 
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Dynamical time: 

Relaxation time due to artificial two-body effects: 



(Binney and Tremaine 1987), where a is the one-dimensional velocity dispersion and (In A) 



9 



is the Coulomb logarithm associated with gravitational interaction, which we estimate by 

(lnA)^ = ln(r/e). (6) 



Cooling time: 



_ 5 nkT 

^cool = 



2 Af-ooi 

= 2.1 X 10''n"^r^/2yr, (7) 
where k is Boltzmann's constant and n and T in the second line are in cgs units. 
• Hubble time (3/2 times the age of the universe): 

tH = H^Hl + z)-^/\ (8) 

Figure |^ shows these timescales as functions of radius at various redshifts. Also shown is the 
heat conduction time tcond) which will be considered in § ^. From Figure |^ several points become 
clear: 

(i) All of the relevant timescales are always larger than (At)sys in any part of the cluster. This 
assures that the time step adopted in our simulation is short enough to accurately follow the 
evolution, even where cooling is taking place. 

(ii) The relaxation time treiax is always much greater than the Hubble time tn- This means that 
the mass resolution is high enough that the MR simulation is not being affected by artificial 
two-body relaxation. 

(iii) The cooling time tcooi becomes shorter than t// at 2 ~ 0.7 both at the center of the main 
clump and at the center of the smaller subclump. This implies that cooling is certainly affecting 
the thermal and dynamical evolution of gas. 



4. Discussion 



We have seen that the high-resolution simulated clusters do not develop a core in the gas 
distribution as is observed. Rather, the gas density increases steeply toward the center. In this 



- 12 - 




Fig. 4. — Dynamical time (long-dashed line), two-body relaxation time (short-dashed line), cooling 
time (solid line) and conduction time (dot-dashed line) in the cluster in run MR at (a) z = 1.4, 
(b) z = 0.7, (c) z = 0.3, and (d) z = 0. The horizontal thin solid line is the Hubble time, and the 
horizontal thin dashed line is the system time step (Ai)sys. The abscissa is the physical radius. 
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section we discuss what is responsible for the discrepancy between the simulated and observed 
clusters. 

To begin with, let us consider numerical resolution. In fact it is uncertain whether we have yet 
achieved sufficient mass and spatial resolution to accurately follow the evolution of the gas density 
profile with bremsstrahlung cooling included. However, even if our resolution is insufficient, this is 
probably not the cause of the discrepancy. If we had better resolution, we would be able to follow 
density inhomogeneities on small scales more accurately. Since the cooling rate is proportional to 
the square of the gas density, this would have resulted in more efficient cooling, and hence an even 
steeper gas density profile toward the center. 

In the remainder of the paper we discuss the effects of (1) changing the cosmological 
parameters, (2) heat conduction, and (3) heating of gas due to supernovae (SNe). 



4.1. Cosmological parameters 

There are now a variety of observations that constrain the cosmological parameters, such as 
h, ^0) ^bi i^S) s-nd 0^, where ilj, is the density of an additional hot dark matter component 
in units of the critical density. It seems reasonable to expect that changing these parameters 
would not largely alter the cluster density profile if the parameters were constrained to satisfy 
observational requirements. 

Since hot dark matter does not cluster on small scales, adding a hot dark matter component 
will tend to make the central density profile shallower. This may push the simulation result closer 
to those of the observations. However, too much of a hot dark matter component results in failure 
to form objects at a high enough redshift, and at the level of currently considered viable models, 
ill, = 0.2 ( Klypin et al. 1995|) , the effects on cluster properties are not large. 



4.2. Heat conduction 

If heat conduction takes place in a timescale shorter than cooling time, it will tend to stabilize 
the gas, perhaps resulting in nearly isothermal gas with a smooth density profile near the center. 
We estimate the heat conduction time icond in the simulated cluster by 

_ 5 nkT 

^^"""^ " 2|V-(kVT)|- ^ ^ 

Here 

K = 1.8 X 10"^r^/^(lnA)e"^ ergss^^cm^^K^^ (10) 
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is the heat conductivity, where the Coulomb logarithm (InA)e, associated with electrostatic 
interactions between electrons, is 



(InA)e = 37.8 + In 



VlQSK 



no 



10 



-1/2- 



(11) 



with rig being the electron number density ( Sarazin 1986| ). 



In Figure ^ we plot tcond as a function of radius. Just outside the radius where tcooi begins to 
drop below tn, tcond becomes smaller than both tcooi and tn- This implies that heat conduction 
may be efficient enough to prevent the inner gas from cooling. 

Thus heat conduction may have played a crucial role in the formation of the core in spatial 
distribution of gas in clusters. However, observations suggest that cluster gas is associated with a 
weak magnetic field, which is tangled on a scale ~ lOkpc. If this is the case, the heat conductivity 
may be suppressed by more than 2 orders of magnitude (Krommes, Oberman and Kleva 1983; 
Chandran and Cowley 1998| ), and conduction may not be able to prevent gas from cooling. At 
present, both the physical state of the magnetic field in clusters (if existent) and the extent to 
which it can suppress conduction remain open questions. 



4.3. SN heating 

Heat input to cluster gas from supernova (SN) explosions occurring in cluster galaxies 
is another process that may be responsible for the discrepancy between the simulation and 
observations. SN heating will raise the entropy and may lower the density and raise the 
temperature of the gas so efficiently as to prevent it from cooling. If this is the case, then the gas 
distribution may develop a core similar to the observed ones. 

To assess this effect, we estimate entropy increase due to SN heating, using observed iron 
abundance. Heat that has been ejected from SNe per unit gas mass can be estimated as 



Aq = 2.0 X lO'^x 



(12) 



where x is the iron abundance in units of solar abundance, Ei and Mpe,! are the explosion energy 
and ejected iron mass per Type la SN, En and Mpe,!! are the same for Type II SN, and /i is the 
fraction of iron that has been produced by Type I SNe. We adopt 

El = 1.3 X 10^^ ergs, (13) 

Mpe,! = 0.744 M0, (14) 

Ell = 1 X 10^^ ergs, (15) 

and 

Mpe,!! = 0.038 Mq (16) 



(Nomoto, Thielemann, and Yokoi 1984; Thielemann, Nomoto, and Hashimoto 199^ ). The value 



of Mpe,!! is an average over the progenitor mass weighted with the Salpeter initial mass function. 
The fraction fi can be estimated from abundance ratios, but there is still considerable uncertainty 
( Gibson, Loewenstein, and Mushotzky 199"^ ). We adopt fi = 0.5. We take the AWM7 cluster as 



an example ( |Ezawa et al. 1997[ ) and fit the observed iron abundance as a function of radius r to 
the following profile: 

r \ 2 



X = 0.58 



^ ^ V0.058/i-iMpc/ 



-0.3 



(17) 



If the SN rate is nearly constant up to the present epoch, we can approximate the entropy 
increase due to SNe, (As)gN, fairly well by Aq/T. Figure |5|a shows (As)sn as a function of radius 
r. We express the entropy in units of 3/c/(2/xmp), where fj, is the mean molecular weight and rrip 
is the proton mass. We also plot s — Si, where Sj and s are the entropy at the initial and present 
epochs in the simulation. In the innermost region of the cluster, (As)sn exceeds s — Si. This 
implies that heating due to SNe can significantly alter the thermal evolution of the gas in the 
central region. 

If the SN rate was much larger at z ^ 1 than at z ~ 0, as is suggested by apparent lack of 
evolution in the iron abundance in clusters toward higher redshifts ( [Mushotzky and Loewenstein 



1997 ), then the above estimate of (As)sn is a poor approximation in the central region because 
the temperature has significantly decreased toward the present epoch. In this case, it is better to 
assess the effect at a higher redshift. Figure ^6 shows (As)gN and s — Si at z = 0.7. This time 
(As)sN is smaller than s — Si everywhere. This indicates that SN heating should have little effect 
in cases in which SN activities in the cluster continued at a nearly constant rate up to z ~ 0.7 and 
have essentially stopped after that. Finally, we can consider the case where most of the entropy 
was injected within the first ~ lO^yr. Then the effect should be very significant, because the gas 
temperature was relatively low: we find n ~ 1 x 10~^cm~^ and T ~ 5 x 10^ K at the density peak 
of a protocluster at z = 5, which will end up in the central region of the cluster at z = 0. Hence 
(As)sN is estimated as ~ 80 and far exceeds s — Sj ~ 4. 

Thus SN heating can be an important factor in the formation of the core in the gas 
distribution, but the effect is strongly dependent on, among other things, the detailed star 
formation history inside the cluster. 



5. Conclusions 

Finally let us briefly summarize our conclusions. 

1. Our multiresolution simulation has followed with reasonable accuracy the evolution of gas 
and dark matter in a typical cluster under the influence of bremsstrahlung cooling as well as 
gravity and hydrodynamics. 
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Fig. 5. — Entropy increase (As)sn due to supernovae {solid line), and entropy change s — in the 
simulation {dashed line), (a) z = 0. {h) z = 0.7. 
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2. The high resolution simulation resulted in a gas density profile steeply rising toward the 
center, with consequent very high X-ray luminosity; however, these properties are not 
observed. 

3. Heat conduction and SN heating are among the processes that may account for the 
discrepancy. Had we allowed for their likely importance in the real world we might have 
been able to recover the observed gas density profile. In a future work we will examine their 
effects by directly incorporating them into the simulation. 
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